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ABSTRACT 

Two-dimensional kinematics of the central region of M 83 (NGC 5236) were obtained through 
three-dimensional NIR spectroscopy with Gemini South telescope. The spatial region covered by 
the integral field unit (~ 5" x 13" or ~ 90 x 240 pc), was centered approximately at the center of 
the bulge isophotes and oriented SE-NW. The Pa/3 emission at half arcsecond resolution clearly 
reveals spider-like diagrams around three centers, indicating the presence of extended masses, 
which we describe in terms of Satoh distributions. One of the mass concentrations is identified 
as the optical nucleus (ON), another as the center of the bulge isophotes, similar to the CO 
kinematical center (KG), and the third as a condensation hidden at optical wavelengths (HN), 
coincident with the largest lobe in 10 /im emission. We run numerical simulations that take into 
account ON, KG and HN and four more clusters, representing the star forming arc at the SW of 
the optical nucleus. We show that ON, KG and HN suffer strong evaporation and merge in 10-50 
Myr. The star-forming arc is scattered in less than one orbital period, also falling into the center. 
Simulations also show that tidal-striping boosts the external shell of the condensations to their 
escape velocity. This fact might lead to an overestimation of the mass of the condensations in 
kinematical observations with spatial resolution smaller than the condensations' apparent sizes. 
Additionally the existence of two ILR resonances embracing the chain of HII regions, claimed by 
different authors, might not exist due to the similarity of the masses of the different components 
and the fast dynamical evolution of M 83 central 300 pc. 
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In spite of the peculiar structure of its circum- 
nuclear region, NGG 5236 presents a well-defined 
bulge that can be fitted by de Vaucouleurs' law 
in a n annular region b e tween radii ^ 10" and 
40" (I Jensen et al.ill98ll: ICallais et al.lll99l[ ) cor- 



responding to « 180 X 720 pc. The dramatic in- 
ward behavior begins approximately at a radius 
of 180-190 pc, where the main dust lanes asso- 
ciated with the galaxy bar spir al into a couple 



of (J-K ) rings that, according to lElmegreen et al 



(Il998l) . might coincide with two inner Linblad res- 
onances. These rings embrace an arc of about 
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twenty star-fo rming condens ations ([Gallais et al 
19911 iHeap etal.lll993l : lHarris et al.ll200lh compa^ 



rable to 30 Doradus in the Large Magellanic Gloud 
Heap et al. 199 3), clu stered in four main knots 
Elmegreen et al. 19981 . Figure 4). The center of 
this arc is approximately w 6.7" to the SW of the 
condensation known as the optical or visible nu- 
cleus (ON) of NGG 52 36. ON marks the center 



of the lElmegreen et al.l (|l998f ) inner nuclear ring, 
but not that of the outer nuclear ring, which is 
shifted 2.5" to the SW of the optical nucleus. K- 



1 



band spectroscopy ( Thatte et al. 2000l ) suggests 
a dynamical center coincident with the center of 
symmetry of the K-band external isophotes (KC), 
which is shifted 4" to the SW of th e optical cen- 
ter (see Figure 3 in Mast et al. 20061 and our Fig- 
ure [T]) . This dynamical center is also suggested, 
although with sm aller spatial resolution , by 2-D 
CO spectroscopy ([Sakamoto et al. I l2004h . The re- 
gion of Thatte's (200 0) isophotal study c oincides 
with that analyzed by I Jensen et al. ( 1981 ). and is 
dominated by the bulge light. 

The disk-like appe arance of the CO kin emat- 
ics interior to 300 pc ( Sakamoto et al. 2004h indi- 
cates that the bar perturbed disk survives well in- 
side the galaxy bulge. How deep it extends into 
the nucleus is a question discussed in this paper. 
The presence of t he hidden condensation (HN) 
(jDiaz et al.ll2006al ). which is larger and more mas- 
sive than the optical nucleus, raises doubts about 
the true nature and the role of the optical nucleus. 
Furthermore, we also cast doubts on the existence 
of a double ILR resonance embracing the arc of 
HII regions. 

To better understand the richness of phenom- 
ena in the central region of NGC 5236 we per- 
formed 3D near infrared spectroscopy at subar- 
second res olution usiiig GEM INI-hCIRPASS con- 
figuration (jPfaz et al.ll2006al ). We observed a re- 
gion including the ON, KC and HN. Our study 
is complemented by N-body simulations to test 
the stability of the M 83 central region, including 
the star-forming arc. In Section [2] we present our 
observations, in Section [3] we discuss the central 
kinematics and the three most important conden- 
sations, which present differentiated kinematics. 
Numerical simulations are presented in Section 3] 
and in Section [5] we present our conclusions. 

2. Observations 

We used th e Cambridge IR P anoramic Survey 
Spectrograph (jParrv et al.l l2000l [2004) installed 
on the GEMINI South telescope in March 2003. 
These observations, performed in queue mode, 
were taken with an integral field unit (IFU) sam- 
pling of 0.36" (6.4 pc) with an array with size 
ss 5" X 13". The array of 490 hexagonal doublet 
lenses attached to fibers provides an area filling 
factor of nearly 100%. The IFU was oriented at 
PA = 120° and centered slightly to the NW of the 



kinematical center according to lMast et al. I (|2002l) 
position. In Figure [H a sketch of the detector po- 
sition is shown. It was exposed for 45 minutes 
and covers a spectral range between 1.2 — 1.4 /xm, 
including Pa/3 1.3 /im and [Fell] 1.26 //m, with a 
spectral resolution of approximately 3200. The 
seeing was about 0.5"; therefore, the focal plane 
was slightly sub-sampled by the configuration con- 
straint. 

The data were reduced using IRAF (distributed 
by the NOAO), ADHOC (2D kinematics analy- 
sis software developed by Marseille's Observatory) 
and SAO (spectra processing software developed 
by the Special Astrophysical Observatory, Rus- 
si a). More de t ails on the reduction are presented 
in iDiaz et al. ( 2006bf ) and the genera l techn ique 
is thoroughly discussed bv lDiaz et al. ( 19991 ) and 
.Mast et al. (2002,) . 

3. 2-D kinematics of M 83 central region 
at less than 10 pc spatial resolution 

In Figure [l] we present the composed JHK 2- 
mass image of M83. The inset shows a blow-up 
of the central 400x400pc, where we detach the 
positions of ON, KC and HN. The astrometric 
positions of the condensations are gi ven in Ta 



ble [H (the astrometry is discussed in iDfaz et al 
2006bf ) . In Figure [21 we show the brightness dis- 
tribution at the Pa/3 continuum, the Pa/3 radial 
isovelocity contours, and the Pa/3 velocity disper- 
sion map. The velocity dispersion is shown as the 
line FWHM (Table [IJ. In spite of the complex 
kinematics shown by the 2-D radial velocity map, 
we can distinguish ordered radial velocities in the 
form of spider diagrams around ON, KC and HN. 
As is well known, a spider diagram indicates disk- 
like motions. The sizes of the disks around the 
three condensations, their maximum radial veloc- 
ities along the line of the nodes, and their velocity 
dispersions at their centers of symmetry are pre- 
sented in Table [H 

3.1. The ionized gas kinematics around 
KC 

KC is the most extended of the three disks. 
Its center coincides, up to the margin of error, 
with the center of symmetry of the bulg e isophotes 



([Jensen et al.lll98ll : iThatte et al.ll2000f ). which in 



turn coincides with that of the CO kinematical 
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Fig. 1. — 2-MASS JHK composite image showing two isophotes and a remarkable coincidence of their center 
with the bulge center KC. The blowup shows a GMOS/GEMINI-S co-added image of the central 20x20 
arcsec with the position of KC, ON and HN. North is up and East is to the left. CIRPASS field true scale 
and orientation are outlined in orange 



center (jSakamoto et al.l 120041 ). We recall that the 
GEMINI spatial resolution is at least twenty times 
higher than that of Sakamoto's CO observation. 

The major axis of the spider diagram is at 
P.A.= 120°. From the spider diagram major- 
to-minor axis ratio, we derive an inclination of 
68° ± 8° for the gaseous disk. We have to be cau- 
tious at this point because the minor axis veloc- 
ity field seems to b e perturbed. The sani e situ- 
ation was faced by I Sakamoto et al. ( 2004 ) , who 
also quote two highly discrepant possible values 
for the inclination of the 300 pc CO disk. These 
authors finally chose the inclination angle of the 



large scale galactic disk based on kinematical rea- 
sons. In our case, the best two-dimensional veloc- 
ity model fitting to the observed velocity field is 
attained for an inclination angle of 25° ± 8° (see 
Section [3.4p . which coincides, up to the margin of 
error, with the small inclination of the large scale 
M83 disk. A small angle also better fits the inner- 
most scale rotation curves to Sakamoto's 300 pc 
disk one (see Fig. 3). The inflow of matter along 
the bar into the central region might be mimicking 
the disks' high inclination. The Pa/3 disk around 
the KC and Sakamoto's CO disk present different 
orientations (see Figure [3]), probably pointing to a 
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Fig. 2. — Panels represent, from left to right, the brightness distribution at the Pa/3 continuum and the Pa/3 
velocity dispersion and isovelocity maps. The Pa/3 continuum map shows the position of ON. The isovelocity 
map shows with crosses the positions of KC in red and HN in yellow. The dashed circles around KG and 
HN represent the astronietric uncertainty in their position. 



phenomenon of tran sition between X2- and a; i -like 
orbits, as quoted bv I Sakamoto et aL ( 20041 ). or to 
a strong per turbation o f the ce ntral structure, as 
discussed by iDiaz et aL 1 (|2006al) . 

The spider diagram (see Figure [2]) can be traced 
up to Umax ~ 65 pc but it is strongly perturbed to 
the E at the receding extreme, and to the W at the 
approaching extreme. The KC rotation curve was 
obtained from approximately 60 radial velocity 
points distributed along the whole disk. From this 
dataset we obtained 13 mean rotation velocities 
between 4pc and 65 pc. Each velocity point en- 
tering the mean was weighted by w{a) = \cos{a)\, 
where a is the position angle of the line joining the 
point to the disk center with respect to the line of 
the nodes. The disk around KC can be described 
in terms of a Satoh's-like spheroid with an effec- 
tive radius of Rkc = 38 ± 8 pc and a mass inside 
Rmax of Mkc ~ (60.0 ± 20) X IO^Mq (for a disk 
inclination of 25°). We note that an inclination 
of 68° for the disk around KC would give a mass 
of « (17 ± 2) X lO^ MfTi, similar to that deduced 
spectroscopically bv lThatte et al. ( 2000[ l from the 



^^CO bandhead velocity broadening. 

The velocity dispersion in the non-resolved cen- 
tral part is of the order of 80kms^^ (see mid- 
dle panel in Figure [2]). Assuming that this veloc- 
ity dispersion is due to a central mass concentra- 
tion, namely a putative BH, we determined the 
BH mass upper limit, Mbh, by adding Keplerian 
rotation curves convolved with a 9 pc Gaussian 
to that of the Satoh disk, and constraining the 
resulting Satoh-|-Kepler rotation curve with the 
central velocity dispersion and the errors of the 
measured velocity points. The result of this pro- 
cedure is shown in Figure[5]for Mbh = 0.2, 1.0 and 
2.0 X 10<^Mq. The best fit is quoted in Tabled 

This analysis points to KC as the true nucleus 
of M83. 

3.2. The ionized gas kinematics around 
HN 

Located at R= 7.8" ± 0.7" to the NW of KC 
(Fig. 2), the spider diagram of the ionized gas 
around HN presents a P.A. = 64°. The form of 
the spider diagram (open legs) points to a more 
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Table 1: Condensation data. The 2a uncertainties in the astrometry of ON, KC and HN are 0.15", 0.8" 
and 0.7", respectively. The three last columns show characteristics of the disk around each one of the 
condensations: ^max is the largest extension along the line of the nodes where disk kinematics can be 
reliably distinguished. Vpeafe is the peak velocity along the line of the nodes, corrected for inclination and a 
is the velocity dispersion at the position of the condensations. 





RA 


Dec. 




^ peak 


a 




(J2000) 


(J2000) 


[pc] 


[kms""""] 


[kms"-'-] 


KC 


13''37™0.56'' 


-29° 51' 56.9" 


65 ±5 


68 ±8 


82± 10 


ON 


13''37"0.98'' 


-29° 51' 55.5" 


8±1 


46 ±9 


110 ± 10 


HN 


13''37"0.46'' 


-29° 51' 53.6" 


45 ±8 


58± 11 


96 ±10 




Fig. 4. — Rotation curves of the Kinematical Center (KC), the Optical Nucleus (ON) and the Hidden 
nucleus (HN). They ar e compared in the right panel to the 300 pc CO disk rotation curve, as obtained from 



Sakamoto et all ( 20041) . 



Table 2: Satoh model- fitting parameters. 





KC 


ON 


HN 


PA [°] 


120 


120 


64 


in 


25.0 ±8 


60± 10 


62 ±8 


Re// [pc] 


38.0 ± 4 


8± 1 


33 ±3 


Mk [10«Mq] 


60.0 ± 20 


4± 2 


20 ±7 


Mbh [10«Mq] 


0.2 - 1.0 


< 1.0 


0.2 - 1.0 



concentrated and more homogeneous distribution 
than that of the KC. It is probably a cannib alized 



dwarf s atellite of M83, as discussed bv Diaz et al 
(I20fl6bl ). To derive the rotation curve as for KC, 
we used a distribution of more than 50 points 
along the whole disk, which provided 13 mean ro- 
tation velocities up to R= 50 pc. In Tables [1] and 
[2] we present the kinematical data for HN and the 



derived Satoh disk parameters. The mass of HN 
was also derived from the rotation curve, as well 
as that of the putative Black Hole at its center 
(see Table [5]), through the velocity dispersion, fol- 
lowing the procedure previously outlined for KC. 

3.3. The ionized gas kinematics around 
ON 

ON is traditionally referred to as the off 



center ed optical nucleus of M83. I Sakamoto et al 



(|2004[ ) proposed that it could be a cannibalized 
satellite. 



10. 



From 



From photometric data. lElmegreeii et al.l (119981) 
derived a mass of 4 X 10^ Mr and iThatte et al 
(l200fll ) of 2.5 X 10^ Mp 



12 



CO 



ban d-head 
( 200nh 

derived a kinematical mass of 1.3x10^ Mm within 



2.293 /im velocity dispersion. iThatte et al 



5 




R[pc] 



Fig. 5. — The rotation curves of (a) KC and (b) HN were fitted by Satoh's models (red line). The dashed 
lines (green, blue and magenta in the electronic edition) were obtained by adding point-like sources of 
Mbh = 0.2, 1.0 and 2.0 x 10^ M©, convolved with a 9 pc Gaussian. Mbh = 0.2 x 10^ Mq is more compatible 
in both cases. 



5.4 pc, assuming that the optical nucleus is viri- 
alized. The mass of ON derived photometrically 
disagrees by a factor of 4 to 5 from that derived 
spectroscopically, probably indicating a system 
with a non-virialized periphery, as discussed in 
Section |4] in connection with the many-body in- 
teraction. The disk around ON is the smallest of 
the three disks discussed in this paper. The ON 
kinematical center is slightly shifted to the SW of 
the ON continuum peak at 1.3 /im, as can be seen 
in Figure [21 The velocity gradient across ON is 
determined from data near the spatial resolution 
limit. We are probably dealing with problems of 
oversampling. The rotation curve with mean er- 
rors is shown in Figure [4] The mass of ON was 
also derived from the rotation curve, as well as 
that of the putative Black Hole at its center (see 
Table [2]), using the velocity dispersion and the 
procedure previously outlined for KC. Our deter- 
mination of the mass agrees very well with that 
obtained from photometric data. 

3.4. Fitting the whole central radial veloc- 
ity field of the gas 

The data in Tables [T] and [2] allow us to gener- 
ate a composed radial velocity map model (Fig- 
ure [6^), to be compared with the observed one. 
The fitting of KC was the most critical due to the 
difficulty in the determination of its disk inclina- 



tion angle ()3.ip . Several fits were made, keeping 
constant the inclination angles of the disks around 
ON and HN and allowing the one associated with 
KC to vary from 25° to 60°. We verified that the 
disk-like residual around KC decreases dramati- 
cally when approaching smaller inclinations. Fi- 
nally, the smallest residuals are obtained for an 
inclination of 25° ± 5°. The final map is shown 
in Figure [5Jd. This inclination is in good agree- 
ment with the inclination of the large scale M83 
disk and that o f the inner 300 pc, d erived from 
CO kinematics ( Sakamoto et al. as can be 



seen in Figure |4l This again points to a continuity 
of the galactic disk kinematics from a kiloparsec 
scale to the smaller scales studied in this paper, 
although variations in the position angles of the 
lines of nodes are observed at 300 pc and at scales 
smaller that 100 pc as well (see Figure [3|). 

Figure [Hb shows that around KC, ON and HN 
the residuals are mainly at the noise level, indi- 
cating that the interpretation of the observed ve- 
locity map as resulting from the superposition of 
three disks is, from a qualitative point of view, 
essentially correct. Nevertheless, deviations from 
rotation points to a more complex situation and 
indicates that the mass estimations based on ro- 
tation curves are barely aproximate. As an ex- 
ample we detect anomalous blueshifted kinematics 
at the East extreme (upper right) above ON, and 
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Fig. 3. — Central part of the I Sakamoto at alj 



(|2004l ) CO (2-1) velocity map (brown contours) 
superimposed on the Pa/3 isovelocity maps. 



redshifted kinematics to the West of HN. These 
regions coincide with the high velocity dispersion 
zones in Figure O These might indicate regions of 
inflow of gas along the bar falling into the central 
region. 

The fit of the rotation curve, from the scale of 
tens to hundreds of parsecs (Fig. 3), does not seem 
to confirm the existen ce of d ouble ILR, as claimed 
( 1998h and I Sakamoto et al 



bv lElmegreen et al 



(120041) 



4. N-body simulations 

In order to understand the dynamical evolu- 
tion of the M83 central region, we performed N- 
body simulations of a system composed of KC, 
HN, ON and four condensations representing the 
star- forming arc. Hernquist's mode ls were used for 
KC, HN and ON (lHernauistJll993h . We included 
stellar and gaseous components. The models were 
constrained by the bodies' rotation curves. Fig- 
ure [7] shows the model rotation curves adjusted in 
order to fit the observational data. The model pa- 
rameters are quoted in Table [3l The four conden- 
sations on the star forming arc were represented by 




Fig. 6. — (a) Velocity map model obtained by 
superposing the disks around KC, HN and KC. 
(b) Residual velocity map resulting from the sub- 
traction of the modeled velocity map from the ob- 
served one. 



Plummer's models, with 1000 particles each, with 
positions and ma sses corresponding to r egions 2 
to 5 in Figure 8 of lElmegreen et al. where: 
region 2 has M2 = 1.3 x 10^ M©; regions 3 and 4 
have M3 = M4 = 1.4 X 10^ M©; and region 5 has 
M5 = 1.8 X 10'' M©. All Plummer models have a 
core radius of 1.5 pc and a cutoff radius of 15 pc. 

A large set of simulations was run in order to re- 
fine the models and orbital parameters used in the 
final simulation. The orbital parameters were cho- 
sen considering that M83 rotates clockwise. Incli- 
nations and position angles are as in Table 2. Here 
we suppose that clusters 2 to 5 on the star forming 
arc, as well as ON are in circular clockwise (pro- 
grade) orbits in the plane of the exte rnal disk of 
M83. According to (jPfaz et al.ll2006al ). HN is the 
remnant of a captured galaxy that triggered the 
star formation along the arc of HII regions. The 
form of the arc and its position with respect to 
HN lead us to adopt a circular counter-clockwise 
orbit for HN. 

The simulat ions were performed with Gadget2 
( Springe] 2005 ). A total of 91,552 particles was 
used, and the simulations were run for 300 Myr, 
starting at the present time configuration. We 
considered starting at a previous epoch, before 
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Fig. 7. — From left to right: Rotation curves of KC, HN and ON modeled with Hernquist compound models. 
Dashed lines are separate components and continuous line is the model circular speed. Black dots with error 
bars are from CIRPASS data. 



the formation of the arc of HII regions, but the 
complexity of the environment at the M83 central 
region and the number of parameters required to 
qualitatively take into account all possible config- 
urations in the past would be a matter for a paper 
especially devoted to that study. 

The simulation (see Figure [5]) shows that the 
galaxy nucleus (KC), the optical nucleus (ON) and 
the hidden nucleus (HN) would form a single mas- 
sive core in about 16 Myr. All Plummer models 
(clusters 2 to 5) fall into the nucleus in 130 Myr. 

Considering the range of uncertainties in the 
orbit determination, we can state that this mas- 
sive core would finally settle as the new nucleus of 
M 83 in a few tens of Myr, implying a net growth of 
the central galactic mass. Furthermore, the whole 
star formation and nuclei merging event would 
last about the time of a global galactic revolution 
(about 150 Myr at the radius of 5 kpc). 

Simulations also show that tidal-striping of the 
condensations boosts the velocity field of the ex- 
ternal shell to escape velocity, which, if mistaken 
for the dispersion velocity of systems in equilib- 
rium, might lead to an ovcrcstimation of their 
masses. This seems to be the case for the dis- 
agreement betwee n the photome t ric an d kinematic 
masses derived by IrhatteeFaD (|2000[ ) for ON. 

The circumnuclear HII regions arc is far from 
being a stable system. In fact, it will spread out 
in an orbital time and be swallowed by the new re- 
furbished nucleus in a period slightly longer than 
the time of merging of ON, KC and HN. There- 
fore, the claimed existence of two inner Linblad 
resonances is difficult sustain. 



5. Conclusion 

We performed CIRPASS 2D spectroscopy at 
Pa/3 and its continuum with the GEMINI-S tele- 
scope. Pa/3 spectroscopy allows for penetration of 
the dust in the direction of KC and HN, and for 
determination of the ionized gas kinematics in the 
central « 5" x 13" embracing ON, KC and HN. 
The kinematics can mainly be explained in terms 
of three disks around each one of the condensa- 
tions, KC, HN and ON. Perturbations are visible 
when the three-disk model is subtracted from the 
observations. The line of the nodes of a tens of 
parsecs-scale disk around KC is rotated with re- 
spect to that of the hundreds of parsecs-scale CO 
one, which in turn is rotated with regard to the 
line of the nodes of the outer galactic disk. They 
are most probably coplanar up to the margin of 
error. Nevertheless, a relative warp of the galactic 
disk inwards of up to 40° can not be ruled out. 

Besides the mass of the condensations KC, HN 
and ON, we inferred upper mass limits for the pu- 
tative BH that could be associated to them. 

Our simulation shows that the three nuclei, KC, 
HN and ON, will merge in a few tens of Myr. The 
disruption of the circumnuclear HII regions arc in 
an orbital time shows that even if ILR do exist, 
this is not enough to ensure the stability of the 
HII regions' orbits. 

We are witnessing a dramatic change in the cen- 
tral region of M83, the nearest galaxy with strong 
nuclear star formation, which may in a few Myrs 
change the mass at its kinematical center by a fac- 
tor of 2, as well as the global aspect at a hundred 
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Table 3: KG, HN and ON model parameters. 





KC model 


HN model 


ON model 


I'l LllliUCl \JL JJCli LHw/iCO ill UioiV 


16384 


8192 


4096 


J— 'iDl^ lliCtOD 


15.0 


13.0 


3.5 


Disk radial scale length 


17.0 


13.0 


4.5 


Disk vertical scale ttiickness 


1.7 


0.8 


0.4 




42.5 


30 


8.0 


Toomrp O at R t 


1.5 


1.5 


1.5 


Number of particles in gas disk 


iDoo4 


Q1 QO 

oiyz 


4uyo 


Gas disk mass 


1.5 


1.5 


0.35 


Gas disk radial scale length 


17.0 


13.0 


4.5 


Gas disk vertical scale thickness 


1.7 


0.8 


0.4 


Toomre Q at Rref 


1.5 


1.5 


1.5 


Number of particles in bulge 


512 


512 


512 


Bulge mass 


1.3 


0.5 


0.3 


Bulge radial scale length 


3.4 


4.0 


1.0 


Number of particles in spherical component 


16384 


8192 


4096 


Halo mass 


150.0 


40.0 


7.0 


Halo cutoff radius 


200.0 


80.0 


45.0 


Halo core radius 


17.0 


25 


4.5 



Note: Simulations were done in a system of units with G=l. Model units scale to physical ones such that: 
length unit is Ipc, velocity unit is 65.58 kms~^, mass unit is 1 x lO^M© and time unit is 14909. 92yr. 



parsec scale. Although HN was interpreted as a of Technology, funded by the NASA and the NSF. 

captured dwarf satellite, the metamorphosis of the 
M 83 central region is independent of the origin of 
HN. 
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Fig. 8. — Evolution of the M83 central region model. ON and the Plummer models (2 to 5) are in prograde 
(clockwise) orbits around KG, while HN is in retrograde orbit. KG, HN and ON merge in 16 Myr. The 
last Plummer sphere (cluster 4) coalesces with the central body in 130 Myr. Gas particles are blue, stellar 
particles red. Spherical components of KG, HN and ON models are not shown for clarity. 
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